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Inelastic Electron Tuneling Spectroscopy (IETS) is a powerful spectroscopy that allows one to 
investigate the nature of local excitations and energy transfer in the system of interest. We study 
IETS for Topological Insulators (TI) and investigate the role of inelastic scattering on the Dirac 
node states on the surface of TIs. Local inelastic scattering is shown to significantly modify the 
Dirac node spectrum. In the weak coupling limit, peaks and steps are induced in second derivative 
d 2 I/dV 2 . In the strong coupling limit, the local negative U centers are formed at impurity sites, and 
the Dirac cone structure is fully destroyed locally. At intermediate coupling resonance peaks emerge. 
We map out the evolution of the resonance peaks from weak to strong coupling, which interpolate 
nicely between the two limits. There is a sudden qualitative change of behavior at intermediate 
coupling, indicating the possible existence of a local quantum phase transition. We also find that 
even for a simple local phonon mode the inherent coupling of spin and orbital degrees in TI leads 
to the spin polarized texture in inelastic Friedel oscillations induced by local mode. 



Introduction: Topological insulators (TIs) are a new 
state of quantum matter, which is insulating in the bulk 
but possesses metallic surface states |H2j . One of the 
defining properties of TIs is their stability to disorder. 
There exists a topological invariant associated with the 
bulk band structure that protects TIs against any time- 
reversal invarant perturbations. Much attention has been 
focusing on impurity scattering on the surface of TIs, and 
a number of novel effects have been revealed, e.g. sup- 
pression of backscattering [3] , absence of Anderson local- 
ization [415] . weak antilocalization [BHS], new quasipar- 
ticle interference patterns [TUIITB] , and impurity-induced 
resonance states [T6HT9] . However the scattering pro- 
cesses considered so far are mostly elastic scatterings. 
Much less effort has been devoted to the study of inelas- 
tic scattering processes that involve not only momentum 
but also energy transfer. Inelastic scattering processes 
are ultimately related to strong correlation effects, which 
is of particular theoretical interest in the study of topo- 
logical matter. 

In this paper, we consider the inelastic scattering of 
Dirac fermions off a single localized vibrational mode on 
the surface of TIs. The local mode is an impurity with 
internal degrees of freedom, similar to a magnetic impu- 
rity. The crucial difference with a magnetic impurity is 
that here time reversal symmetry is preserved. A related 
model of a Fermi gas interacting strongly with a localized 
phonon was proposed by Yu and Anderson [20] , where it 
was found that the Fermi gas softens the phonon fre- 
quency and generates an effective double well potential 
for phonons at sufficiently strong coupling. This model 
has been heavily investigated subsequently (see e.g. [2TT 

W)- 

The local electronic structure is also substantially mod- 
ified by the presence of the local mode, which can 
be detected by inelastic electron tunneling spectroscopy 
(IETS) with a scanning tunneling microscope (STM) 
[25126] . The pre-STM IETS using metal-insulator- metal 



tunneling junctions measures the vibrational spectrum 
of the dynamical defect on the insulator, with the peak 
in d 2 V/dI 2 corresponding to the vibrational frequency 
[37H3D]. IETS for s-wave superconductors was consid- 
ered in [31], and for d-wave superconductors in [32H34] . 
where it was found that at weak coupling there is a kink 
in the dl/dV curve |32j . and at stronger coupling reso- 
nance peaks emerge [33] . The interference of incoming 
and outgoing waves in inelastic scattering can also pro- 
duce standing wave patterns, i.e. inelastic Friedel oscil- 
lations [35136] , 

In this paper, we calculate the local density of states 
(LDOS) for Dirac fermions in the presence of a local 
mode. The main results are as follows: (i) In the weak 
coupling limit, by using second-order perturbation the- 
ory, we find that at the impurity site, there is a kink and 
a logarithmic singularity at the local mode frequency u>q 
in the fermion LDOS. Correspondingly in the frequency 
derivative of LDOS, i.e. d 2 I/dV 2 , there is a step and 
a peak at ujq. Away from the impurity site, inelastic 
Friedel oscillations are shown to be present in d 2 I/dV 2 . 
Similar features are found for spin polarized STM. (ii) 
In the strong coupling limit, where we have a single site 
problem, the LDOS becomes a series of delta functions, 
forming a single band at negative frequency. The density 
of states at the Fermi level is fully depleted. The result 
here is universal in the sense that it is independent of 
the band structure of conduction electrons, (iii) By sum- 
ming over the noncrossing diagrams, we also map out 
the evolution of LDOS from weak to strong coupling. At 
intermediate coupling two resonance peaks appear. As 
coupling further increases, the peak at positive frequency 
first becomes higher than the one at negative frequency, 
and then the two peaks change role, with the negative 
frequency peak higher than the positive frequency one. 
At even stronger coupling, the positive frequency peak 
finally dies out. 

Inelastic scattering with a localized mode: The 
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FIG. 1: Change of LDOS and its frequency derivative at the 
impurity site due to the presence of a local vibrational mode 
in the weak coupling limit. Here T = 0, frequency is measured 
in units of meV, with local mode frequency ojo = lOmeV, and 
the cutoff A = 300meV. 
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FIG. 2: IETS maps of d 2 5I{r, V)/dV 2 at the voltages (mV) 
indicated with T = 1 K, for spin unpolarized STM (up) 
and for SP-STM with spin-polarization given by mo = 
3n(l,0,0)/4 (down). 



low energy Hamiltonian describing the massless Dirac 
fermions coupled to the local boson at r = is 

H = ^2e k cl a c ka +uj blb +g cj^cv^Oo + b o), (!) 



ka 



kk' a 



with £fc = VfZ ■ (o - x fe), where Vf is the Fermi velocity, 
c = ("a;) GyiGz) is the vector of Pauli matrices. The bare 
fermion Green's function Go (fe, iu) n ) — l/[iuj n —v f{a x k y — 
<j y k x )] reads in real space and real frequency 
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where Ji and Yi are the Bessel function of the first and 
second kind. The local bosonic mode has propagator 
T>(k,ioj n ) = V(iuj n ) = 2w /[(iw„) 2 - w 2 ]. 

Multiple scattering of Dirac fermions off the lo- 
cal bosonic mode can be accounted for by the 
T-matrix, which modifies the fermion Green's 
function to G{r,r'\iu) n ) = G (r,r';iuj n ) + 
Go(r, 0; zw ra )T(iw n )Go(0, r'\ iui n ). The T-matrix is 
of the form T(«w) = E(iw n )/[1 — Go(0, 0; £w n )E(iw„)], 
where the local self energy is calculated from the 
one particle irreducible (1PI) diagrams, with result 
E(iw„) = - 5 2 T£ m G(0,0;zw„ - m m )v(in m ). 



Weak coupling limit: When the coupling between 
the fermions and the local mode is weak, perturbation 
theory can be employed. In this case, the change of 
LDOS is small. However, sharp features can be detected 
in the frequency derivative of LDOS. 

For weak scattering, the T-matrix can be approx- 
imated by the one loop self energy, Ti^°'(iui n ) — 
— g 2 TJ2 rn Go(0,0;iw„ — iQ m )T>(iD, m ), which reads at 
zero temperature 



E(w) = j 

v f 
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+in [(cj - uj )Q(uj - ujq) + (u -> -w)]} , (3) 



where A is a cutoff. At high temperatures, the imaginary 
part of the self energy at zero frequency becomes linear 
in temperature, i.e. ImE(w = 0, T) ~ — (2n 2 g 2 / v^)T . 
It is interesting to compare this result with the effect of 
inelastic scattering in normal Fermi liquids, where the 
self energy has T 2 dependence due to Pauli blocking in 
the absence of a finite volume Fermi surface. 

For Dirac fermions coupled to a localized magnetic 
atom via exchange interaction JJ^S • c ka crck>a, the T- 
matrix becomes g 2 T^2Go(k 1 iu> n — iO, m )aiajXij(i^m)j 
where Xij is the local spin susceptibility. We notice that, 
for an isotropic Xij> the effect on the LDOS will be the 
same as that of a local phonon at least for weak coupling. 

The effect of such inelastic scattering processes 
can be detected by STM, where LDOS N(r,u) = 
— (l/7r)TrIm G(r, r; u>) is measured. The change of 
LDOS has two contributions, one is of the form SN^- 1 ' ~ 
Rc [TrGg(r,0;w)] ImE(cj), the other reads 5N ( ^ ~ 
Im [TrGj^r, 0; J)] ReE(w). 

Consider first LDOS at the site of the local vibrational 
mode. Noticing that G o (0, 0; uj) = ^ [^f In M 
for lu > 0, T = 0, one obtains 



5N {2) = -\uj 2 C-2 



5N {1) = Aw 2 Ci(cj- u )Q(u- w ), (4) 
(w - w ) In - — - — r + (w -> -w ) 



|w - w 1 
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where d = 7r[(4/7r 2 ) ln 2 (w/A) - 1], C 2 = (4/tt) ln(w/A), 
and the overall coefficient A cx g 2 . Such a signal can be 
seen most clearly from the frequency derivative of LDOS, 
i.e. d 2 I/dV 2 (see Fig. [l]). Only above the threshold 
of the boson frequency can the local mode be excited, 
scattering over which renders the fermion lifetime finite. 
Thus ImE(w) is nonzero only for w > wo. This gives rise 
to a kink at w = w in SN^\ or equivalently a step in 
duSN^ 1 '. Such a sharp feature also enters ReE(w) via the 
Kramers-Kronig relation, resulting in a peak in d^SN^K 
The effect of inelastic scattering will be present even 
away from the scattering center. Interference between 
incoming and outgoing Dirac fermions produce standing 
waves on the surface of a TI, similar but qualitatively 
different from conventional Friedel oscillations. These 
standing waves can be seen in the oscillations of the sec- 
ond derivative d 2 I(r, V)/dV 2 , in the narrow window of 
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energies near the energy of the bosonic mode luq where 
inelastic scattering occurs. 

From SN(r, lu) ~ Gq(v, 0; lu)Y,(lu), one can see that the 
spatial dependence of the change of LDOS comes from 
the free fermion Green's function Go(r,0;ui). From the 
asymptotic behavior of the Bessel functions, 5N(t,lo) ~ 
(1/r) cos(2cor / v f +60), the period of the standing waves 
is set by the bias and the Fermi velocity. For the on-site 
LDOS, the sharp features, which have their origin from 
S(o;), occur near the boson frequency luq. We plot in 
Fig. [2] the spatial variation of d w 6N(r,ui) at T = 1 K. 
The main change in this variation as one sweeps through 
the inelastic mode is the center peak (below luq) which 
turns into a dip (above wq), 

A salient feature of the surface state of topological in- 
sulators is the strong coupling between spin and orbital 
degrees of freedom, resulting in the locking of spin dy- 
namics to that of charge. Such a feature is also man- 
ifest in inelastic scattering processes. Even scattering 
off a non-magnetic mode produces sharp features in the 
magnetic structure M(r,u) = — Trim G(r, r; u;)cr/(27r), 
which can be detected by spin-polarized STM (SP-STM) 
with a magnetic tip, distinguishing itself completely from 
a normal metal, for which M(r,u>) vanishes identically 
in the presence of non-magnetic impurities. 

The second order change of M{r,uS) due to the pres- 
ence of a local mode is 6M(r,cj) ~ lmg 'Sg±[f x z], 
where Gq(t,w) = go(r,u))a° +g±(r,u))[rxz]-cr. One can 
see immediately that SM Z = 0, i.e. the induced magnetic 
structure is in-plane only. We plot in Fig. [2] the spatial 
variation of du,SM(r, lu) at different frequencies, with the 
magnetic direction of the tip parallel to the surface. One 
can clearly see anisotropic inelastic Friedel oscillations 
near the boson frequency luq. 

Besides local probes, e.g. STM, the effect of inelas- 
tic scattering on the surface of TI can also be mea- 
sured by more averaging probes, e.g. ARPES, trans- 
port. To illustrate this, we consider a finite but low 
density of localized vibrational modes on the surface of 
a TI. Up to first order in the impurity density n- lmp , 
the impurity-averaged fermion Green's function reads 



G(k,w) = G (k,u) l-S(w)G (fc,w) 



, where the self 



energy in the Born approximation is = n{T{iS) ~ 

niE(w). The resulting spectral function A a {k<jS) — 
-(l/TT)lmG aa (k,Lu) at T = is plotted in Fig. [| For 
|w| < luo, since ImS(cj) = 0, the spectral function is not 
smeared by inelastic processes.^ Only the dispersion is 
slightly modified, with A(k,u) = (l/2)[5(w - ReE(cj) + 
vpk) + 8{lu — Re£(u;) — vpk)] 1 retaining the Dirac cone 
structure. For \lu\ > lu , inelastic scattering gives rise to 
finite ImS(w), and the spectral function is broadened. 
There is a qualitative change at \lu\ = luq. Although 
dressed by thermal fluctuations, one can see from Fig. [3] 
that such a pattern is still visible at relatively high tem- 
peratures. 

Strong coupling limit: We consider here the strong 
coupling limit, where two fermions with opposite spins 
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FIG. 3: Impurity-averaged spectral function at T = 0, 23K 
for k y = 0, with lu = 10meV,A = 300meV, TTg 2 ni/v 2 f = 0.1. 



are trapped on the impurity site. In this case, one only 
needs to consider the single site problem with coupled 
fermion and boson. This is a standard textbook problem 
(see e.g. [37]). However the problem here has a crucial 
difference with the independent boson model solved in 
[3 7) . Here since the effective onsite energy for fermions is 
negative, in the ground state the fermionic state is doubly 
occupied. 

As g — > 00, the model reduces to a single site problem, 
H = J2a t<ycic a +uj tfb+gY Ja clc a {tf+b), which can be 
exactly diagonalized by making the transformation H = 
e s He~ s , with s = ^ c^c CT (6' t — b), and the resulting 
Hamiltonian reads 



H = LUntfb 



E 



LUQ 



— c}c T c} Ci . (6) 



The last term represents an attractive Hubbard interac- 
tion generated by electron-phonon coupling. The local 
mode thus forms a negative-U center. 

With e a = 0, in the ground state, two electrons are 
trapped at the impurity site, while the phonon occupa- 
tion is zero. We also notice that, in the strong coupling 
limit, the effective on-site energy e a —g 2 /LUo is always neg- 
ative, independent of the band structure. The real-time 
electron Green's function G aa '(t) = — i{Tc a {t)c^, (0)) 
can be written in terms of the transformed quantities, 

G aa >{t) = t9{-t)e^TT^e-^cle im c n e' tm ). Here 



c,rX and cL 



and the operator X 



exp 



Since the transformed Hamiltonian is 

a summation of the electron part and the phonon part, 
H = H c + Hb, the Green's function can be factorized as 



e- p6 >Xl(0)X(tj\ , (7) 



The phonon part gives contribution e *W, with Q(t) — 

^2 (l — e 4 " *) at zero temperature. The electronic part 

leads to exp \—it (e CT — ig 2 /wo)] ■ The resulting spec- 
tral function is a series of delta functions forming a single 
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FIG. 4: LDOS at the impurity site in the strong coupling 
limit, with g — 22 (left) and g — 8 (right). Here ujq = lOmeV. 



band, 



nl 



n=0 



UJQ 



- nuiQ 



(8) 



with A = g 2 /uj 2 . 

Evolution from weak to strong coupling: Hav- 
ing considered two limiting cases with weak and strong 
coupling, now we proceed to study generic coupling 
strength, where further approximations are required, and 
the LDOS can only be calculated numerically. We will 
use the self consistent Born approximation, where the 
T-matrix obeys the integral equation, 



r- 1 = 



y(0) , q- 
n ' / j ^"nm ' m 



_ G (o) 



(9) 



where G%> = G o (0,0;^„) = (-*w„/47re^)ln(A 2 /^), 
and /C nm = -g 2 TT>(iui n - iv m ) [G o (0, 0; iv m )\ . 

The numerical results for T = are shown in Fig ^ for 
various coupling strengths. As we increase coupling, the 
LDOS starts to deviate from the Dirac cone strucutre at 
uj > u}q (Fig [5^,), with the spectral weight accumulating 
at intermediate frequencies. Two resonance peaks form 
as coupling further increases (FigJSJs). At still stronger 
coupling (Fig [BJ;), the two resonance peaks move towards 
lower frequency, with the peak height increasing sharply, 
the peak width becoming much narrower, and the asym- 
metry in two peaks more pronounced. 

We notice that in Fig (J5^,,b,c), the linear dispersion is 
retained near 10 = 0. The LDOS is modified only at larger 
frequencies. However, qualitative changes occur at even 
stronger coupling. While from Fig ^jp to Fig [5];, more 
and more spectral weight moves into the peak at positive 
frequency, in Fig [5ji, the peak at negative frequency be- 
comes the dominant one. With increasing coupling, the 
peak at positive frequency is further suppressed (Fig[5^), 
and finally becomes invisible (Fig [5]:), as expected from 
the above result in the strong coupling limit. The LDOS 
at negative frequency also gets depleted near oj = 0, with 
the Dirac cone structure fully destroyed (Fig[5jl,e,f). 

While the changes from Fig [5^, to Fig [5J;, and from 
Fig[5ji to Fig[Hf are continuous, there is a sudden quali- 
tative change from Fig[5j: to Fig[5ji, indicating the possi- 
ble existence of a phase transition 42 . More sophisticated 
methods, e.g. numerical renormalization group (NRG), 
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FIG. 5: (Color online) Evolution of LDOS from weak to 
strong coupling. Here cjo = lOmeV, A = 300meV,w/ = 1, 
and g = 0.003, 0.01, 0.03, 0.1, 0.3, 1. 



are needed to settle this issue. The model considered 
here can be regarded as a generalization of the exten- 
sively studied pseudogap Kondo model (PKM) (see [55] 
and references therein), with the defects having discrete 
internal degrees of freedom generalized to defects hav- 
ing continuous internal degrees of freedom. In PKM, 
where the conduction electrons have a power law DOS 
p(e) ~ |e| r , weak coupling renormalization group calcu- 
lations found an unstable RG fixed point, indicating a 
phase transition between weak and strong coupling |39) . 
However this method is only valid for small r, and non- 
perturbative NRG studies showed that the phase transi- 
tion disappears for r > 1/2 in the particle-hole symmetric 
case [4"0] . 

Conclusion: In conclusion, we have studied the 
inelastic scattering of Dirac fermions with a single im- 
purity that has continuous internal degrees of freedom 
and does not break time reserval symmetry. For weak 
coupling, sharp features near the local mode frequency 
(step and peak) can be detected in d 2 I/dV 2 using IETS- 
STM. At intermediate coupling, resonance peaks appear 
in dl/dV curve. The original Dirac cone structure is 
fully destroyed locally at strong coupling. There is possi- 
bly a zero temperature phase transition, and thus a local 
quantum critical state, at intermediate coupling, which 
we leave for future work. The presence of a finite density 
of such local modes will substantially change the trans- 
port properties of TIs. It would be interesting to explore 
the dephasing effect of such modes, to study how weak 
antilocalization is influenced by inelastic scattering. 
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